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INTRODUCTION 

Plant  perturbations  can  destabilize  a  nominally  stable  system. 
Therefore,  the  term  robustly  stable  refers  to  the  extent  to  which  a 
model  of  the  open-loop  system  may  be  changed  from  the  nominal 
design  without  destabilizing  the  overall  closed-loop  feedback  system. 
What  is  really  desired  is  a  robustness  analysis  that  will  apply  to 
simultaneous  independent  and/or  dependent,  not  necessarily  small, 
perturbations.  There  is  a  certain  range,  due  to  neglected 
nonlinearities  and  unmodeled  system  dynamics,  where  the  model 
and  system  may  behave  in  grossly  different  ways.  Unfortunately, 
this  range  is  implicit  in  the  technology  that  the  model  approximate?, 
so  any  general  theory  must  encompass  a  variety  of  perturbations. 
The  importance  of  obtaining  robustly  stable  feedback  control  has 
long  been  recognized  by  designers.  In  fact  for  standard  single-input, 
single-output  (SISO)  systems,  the  gain  and  phase  margins  are 
basically  measures  of  how  close  the  Nyquist  plots  come  to  encircling 
the  point  -1  [Referencel,  p.  256].  Even  for  SISO  systems,  gain  and 
phase  margins  do  not  exhaust  all  the  robustness  issues.  General 
robustness  measures  are  a  significant  theoretical  challenge  for 
multiple-input,  multiple-output  (MIMO)  control  systems. 

Recently,  a  promising  analysis  tool  called  structured  singular 
values  was  proposed  by  Doyle  {Reference  2]  and  Doyle,  et  al 
[Reference  3]  that  can  be  used  to  nonconservatively  calculate 
stability  margins  in  feedback  loops  with  real  or  complex  variations. 
A  key  observaton  in  Doyle’s  structured  singular-value  theory  is  that 
a  linear  inter-connection  of  inputs,  outputs,  transfer  function, 
parameter  variations,  and  perturbations  can  be  rearranged  to  isolate 
all  perturbations  in  a  feedback  block  similar  to  Figure  1. 
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FIGURE  1.  Perturbation  Feedback. 


The  transfer-function  matrix  M  represents  the  nominal  plant  and 

the  block-diagonal  matrix  A  represents  the  perturbations.  The  block- 
diagonal  elements  are  bounded,  but  otherwise  unknown  real  or 
complex  variables.  The  papers  by  deGaston  [References  4  and  5] 
include  a  highly  informative  example  illustrating  this  technique. 

The  state-space  perspective  introduced  by  Doyle  and  Packard 

[Reference  6]  includes  an  extension  of  structured  singular  values  to 
time-varying/nonlinear  systems.  The  goal  of  this  report  is  to 
introduce  a  state-space  model  for  real-  and  complex-additive 
perturbation  errors  that  is  algebraically  and  dynamically  equivalent 
to  Doyle's  structured  singular  value.  The  frequency-domain  and 
state-space  models  can  be  used  interchangeably  because  they  are 
equivalent.  All  of  the  proofs  and  equations  in  this  report  are 

elementary,  but  the  exposition  does  provide  a  record  of  the  exact 
link  between  the  two  models.  The  fact  that  the  models  are 
equivalent  is  not  surprising.  However,  the  exact  equation  and  form 
of  the  models  that  make  them  equivalent  are  perhaps  not  immediate. 
Moreover,  the  state-space  equations  add  a  degree  of  flexibility  to 
robustness  analysis  that  allows  the  designer  to  use  the  most 
convenient  approach.  A  second-order  oscillator  model  with  a 
product  term,  an  alternative  and  simple  proof  of  the  maximum 
principle,  and  a  Monte  Carlo  search  algorithm  are  all  included  as 
examples  that  support  the  utility  and  flexibility  of  state-space 
models.  These  models  further  demonstrate  that  structured  singular 
values  are  a  convenient  and  natural  way  to  model  additive  errors  in 
linear  systems. 

The  feedback-perturbation  model  is  based  on  the  observation 

that  additive  errors  can  be  represented  as  a  matrix  product.  The 
product  can  be  interpreted  as  the  product  of  a  control  and 
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observation  matrix  with  A  sandwiched  in  between.  In  state-space 
terminology  the  problem  can  be  described  as  a  spectral  assignment 
problem  using  output  feedback,  when  the  variable-feedback  matrix 
A  is  constrained  to  the  block-diagonal  elements  represented  in  Figure 
1.  Combining  this  theory  with  some  well-known  determinant 
identities,  the  output-feedback  matrix  formulation  is  equivalent  to 
Doyle's  feedback  block  diagram  (Figure  1).  Using  the  Sherman- 
Woodbury-Morrison  formula  (Reference  7),  the  output  additive- 
error  model  leads  to  a  perturbed  dynamical  system  that  can  always 
be  decomposed  into  a  parallel  interconnection  of  a  nominal  and  a 
perturbed  system. 

Next,  a  link  between  root  locus  and  structured  singular  values 
is  established.  When  this  link  is  applied  to  a  third-order  polynomial, 
the  robustness  estimates  predicted  by  root  locus  and  structured 
singular  values  can  be  significantly  different.  A  second-order 
oscillator  model  with  product  terms  shows  that  Doyle's  block  diagram 
(Figure  1)  does  not  include  all  additive-error  models  even  for  linear 
systems.  However,  the  state-space  theory  in  this  report  does  admit  a 
larger  class  of  perturbation  models  that  can  accommodate  product 
terms. 

To  apply  this  theory,  a  Monte  Carlo  search  algorithm  was 
combined  with  some  graphical  convergence  strategies  by  Hewer, 
Klabunde,  and  Kenney  [Reference  8]  to  compute  p  for  third-order 
systems.  The  Monte  Carlo  algorithm  was  also  used  to  compute  the 
structured  singular  values  for  the  autopilot  example  [Reference  9] 
first  studied  by  de  Gaston  and  Safonov  with  general  numerical 
agreement.  For  the  de  Gaston  example,  the  Monte  Carlo  algorithm 
yields  sharper  bounds  than  the  Fan-Tits  algorithm  [Reference  10] 
and  Doyle's  estimates,  which  are  consistent  with  the  well-known  fact 
that  their  estimates  are  valid  for  both  complex  and  real 
perturbations. 

While  both  real  and  complex  perturbations  can  be  used  to 
model  additive  errors,  the  role  of  complex  perturbations  can  be 
extended  to  include  dynamical  perturbations.  These  perturbations, 
which  must  be  dimensionally  compatible  with  the  block-diagonal 
structure  outlined  in  Figure  1,  represent  a  family  of  neglected 
dynamics.  Usually  the  family  is  not  included  in  the  initial  design 
because  it  was  considered  an  inessential  component.  However,  the 
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impact  on  the  design  robustness  is  essential.  The  family  of  neglected 
dynamical  systems  is  interconnected  to  the  original  system  by 
feedback.  Within  this  framework  the  original  system  is  converted 
into  a  higher-order  dynamical  system.  A  stability  theorem  for  the 
higher-order  dynamical  system  is  established  in  the  Perturbation 
Dynamics  section. 

In  the  section  entitled  New  Upper  Bound  for  p,  a  new  upper 
bound  for  the  structured  singular  value  is  introduced.  Depending  on 
the  matrix  norm,  this  bound  is  determined  by  a  Kronecker  sum  or  by 
the  solution  of  a  Lyapunov  equation. 


STRUCTURED  SINGULAR  VALUE 

In  order  to  introduce  Doyle’s  structured  singular  value  p,  the 
block  structure  A  in  Figure  1  and  related  concepts  are  defined.  The 

set  of  mxn  real  (complex)  matrices  is  denoted  by  Rmxn(cmx',)-  The 
direct  sum  [Reference  11)  of  two  matrices  E  ©  FOEe  C'1*",  F  e  C'n*m)  is 
the  (n  +  m)  x  (n  +  m)  block  matrix 

rE  °i 

Lo  Fj 

with  compatibly  dimensioned  null-matrices  on  the  crossdiagonal. 
The  direct  sum  of  three  or  more  matrices  follows  a  similar  pattern. 

The  Frobenius  norm  of  a  square  matrix  Ee  C'1*"  is  defined  by 

II F  ||  F  =  (trace  E*E)1/2 

(*  is  the  conjugate  transpose)  and  the  spectral  norm  or  2-norm  of  E  is 
the  square  root  of  the  maximum  eigenvalue  of  ElE,  wherelxll 
represents  the  2-norm  of  the  vector  x 

II E  ||  =supremum  ||  Ex  || 
llx||=l 

If  the  distinction  is  unimportant,  then  11  will  represent  either  norm. 
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The  block-diagonal  matrices  A  in  Figure  1  are  members  of  a  set 
A,  which  represents  the  perturbation  structure  that  is  dimensionally 
and  functionally  compatible  with  the  transfer  matrix  M  (s).  Here  s  is 
the  Laplace  transform  complex  variable  with  imaginary  part  /co 
(j  =  and  to  e  R).  The  operator  Res  represents  the  real  part  of  the 
complex  variable  s.  The  elements  of  A  have  block-diagonal  structure 

A  =  {  A  =  ©£1  A;  }  (2.1) 

where  each  ki  x  kx  block  matrix  Aj  is  further  required  to  have  some 
specific  entries  equal  to  zero  while  other  entries  are  free  to  vary  as 
complex  or  real  numbers.  The  set  A  is  closed  in  the  space  of  (ki  +  k2  + 
...  +  km)  x  (ki  +  k2  +  ...  +  km)  dimensional  matrices  using  the  metric 
induced  by  the  2-norm  or  F-norm  topology. 

This  block  structure  [kj,  k2,  ....  km  ]  is  equivalent  to  the 
decomposition  of  the  perturbation  interconnection  into  parallel 
subsystems.  Nonsquare  perturbations  can  be  accommodated  by 
augmenting  the  interconnection  structure  with  rows  and  columns  of 
zeros.  Reducing  the  members  of  A  to  a  uniform  norm  bound  requires 
that  nonsingular  scaling  matrices  be  absorbed  into  the 
interconnection  structure  in  Figure  1.  The  term  block  structure  and 
the  symbol  A  follow  Doyle's  usage. 

Frequently,  the  block  matrices  Aj  with  block  structure  fkj,k2, 
...,  km]  are  required  to  have  an  even  simpler  structure.  Especially 
useful  is  the  form  of  a  scalar  times  the  identity  matrix:  A4  =  Cj  I.  Let 
D  denote  the  set  of  As  where  each  block  matrix  has  this  form 

£>  =  (A  =  ©£,  Cj  I,  Cj€  CorR}  (2.2) 

This  set  will  sometimes  be  denoted  by  D(cj,  ...,  cm).  We  will  assume 
throughout  unless  otherwise  noted  that  the  structural  constraints  on 
the  elements  of  A  are  such  that  D  c:  a. 

For  a  given  transfer  matrix  M ,  define  the  matricial  spectral  set 
with  respect  to  A  by 
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MS  (M,A)  =  (Ae  A:  det  fl  -  M(s)A)  =  0  for  some  s,  Res  £0}  (2.3) 

The  set  MS(M,A)  may  be  empty,  especially  if  the  elements  of  A  are 
required  to  he  real  matrices.  For  a  constant  nxn  matrix  P,  we  use  the 
symbol  *{P)  to  denote  the  set  of  eigenvalues  of  P  and  define  the 
spectral  radius  of  P  by  p  (P)  a  max  { Ixf,  X  e  A(P)}.  Classifying  the  set 
A  and  the  transfer  matrices  M  for  which  MS  (M,A)  is  nonempty  is  an 
example  of  the  following  inverse  multiplicative  eigenvalue  problem. 
Friedland  [Reference  12]  proved  that  given  a  matrix  Ee  C"*"  and  a 
specified  set  of  eigenvalues,  I,  there  always  exists  a  complex  diagonal 
matrix  D  such  that  A  (ED)  =  £,  provided  all  the  principal  minors  of  E 
are  different  from  zero.  However,  a  general  existence  theorem  for 
the  inverse  multiplicative  eigenvalue  problem  for  real  matrices  is 
apparently  unknown.  As  will  soon  be  demonstrated,  Doyle's  inverse 
multiplicative  eigenvalue  problem  is  equivalent  to  the  well-known 
problem  of  arbitrary  pole  assignability  by  static-output  feedback. 
The  recent  paper  by  Reinschke  [Reference  13]  and  associated 
references  summarizes  recent  work  on  the  pole-assignment  problem 
by  output  feedback.  Reinschke's  algebraic  tests,  which  guarantee 
eigenvalue  placement  using  output  feedback,  while  relevant  and 
illuminating,  do  not  completely  characterize  the  nonempty  sets 
MS(M,A)  for  real  As.  Doyle  defined  the  structured  singular  value  of  M 
with  respect  to  A  by 


p(M,  A)  = 


supremum 


Ae  MS(M,  A  )  11 


0 


if  MS  (M,  A)  *  4> 


if  MS  (M,  A)  =  <j>  (2.4) 


This  optimization  problem  will  generally  have  multiple  local  maxima 
that  are  not  global.  Some  exceptions  for  the  block  structure  [kj,  k2, 

kj]  with  the  complex  perturbation  set  A  were  identified  by  Doyle  and 
further  characterized  by  Fan  and  Tits.  Nevertheless,  Doyle  derived 
the  following  bounds  for  p.  for  a  fixed-block  structure  [kj,  k2,  ...»  km] 
and  D  c  a 


p(M)£  p(M,A)  =  h(DMD-1,A)£||m||2 


(2.5) 
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where 


p(M)=  sup  p(M(jto))  and  ||m||2  =  sup  ||  M(/co)  ]|  2 
o)g  R  toe  R 

When  A  is  complex,  then  some  matrix  in  the  set  D  is  always  a 
member  of  MS  (M,  A),  and  p(M)  is  always  a  lower  bound  for  p.  When 
A  is  real  and  the  matricial -spectral  set  is  nonempty,  then  p(M)  is  still 

a  lower  bound  for  p  provided  some  member  of  the  set  D(X,  X . X)  for 

real  X  is  in  MS  (M,  A),  because  the  matricial-spectral  set  over  real 
perturbations  is  always  a  subset  of  the  matricial-spectral  set  for 
complex  perturbations. 


LINEAR  SYSTEM  MODEL  FOR  ADDITIVE  ERRORS 

In  this  section  a  state-space  model  of  additive  perturbations  is 
developed.  First  the  standard  linear-control  system  is  defined,  along 
with  state-space  counterparts  of  the  matricial-spectral  set  and  the 
structured  singular  value  p.  Next,  Theorem  3.3  is  proved,  which 
shows  that  we  may  restrict  our  attention  to  the  imaginary  axis  when 
computing  elements  in  the  matricial-spectral  sets.  The  main  result  of 
this  section  is  Theorem  3.4,  which  establishes  the  equivalence  of  the 
frequency-domain  and  state-space  matricial-spectral  sets,  and  their 
associated  structured  singular  values. 

Consider  the  fallowing  linear  dynamical  system 


x  =  Ax+  IB),  B2] 


- 1 

"Ci 

Ly2. 

_C2_ 

(3.1) 

(3.2) 


The  matrices  A,  Bj,  B2,  Cj,  and  C2  are  nxn,  nxl,  nxk,  pxn,  and  kxn 
dimensional  real-valued  matrices,  respectively.  The  state  vector  x, 
the  input  vectors  rj  and  r2,  and  the  output  vectors  y\  and  y2  are 
compatibly  dimensioned. 


After  taking  the  Laplace  transform,  we  obtain  the  two  port 
with  (p+k)x(l+k)  block  transfer  matrix  T(s) 


9 


NWC  TP  6971 


Yll  =  Hs) 

Ri")  _ 

Gll(s)  Gl 2(s) 

Ri' 

[Y2. 

,R2j  = 

_G21(s)  G22(s). 

.R2. 

The  triples  (A,  Bj,  Cj)  and  (A,  B2,  C2)  are  both  minimal  (i.e., 
controllable  and  observable)  realizations  of  their  respective  transfer 
matrices 


Gij(s)  =  Ci  (si  -  A)  !Bj  i  =  1,2,  j  =  1,2 

The  plant  matrix  A  incorporates  both  open-loop  plant  dynamics  and 
any  feedback  compensation  employed.  We  assume  throughout  this 
paper  that  A  is  stable,  that  is  the  spectral  abscissa  a  (A)  is  negative 

a(A)s  max  {ReX:  X  e  A(A)} 


Perturbations  of  A  of  the  form  A  +  8A  will  represent  additive- 
perturbation  errors  in  the  dynamical  system  (7  1).  We  assume  that 
the  matrix  term  8A  can  be  factored  into  the  matrix  product  6A 
=  B2SAC2  for  some  A  €  A.  Thus,  the  additive  perturbation  matrix  5A 


is  obtained  by  applying  the  output  feedback  vector  r2  =  SA  y2. 

f  m  \ 


block  structure  parameters  [kj,  k2,  ...km] 


k  =  E  k: 
i=l  lJ 


characterize 


The 

the 


additive  perturbation’s  coefficients  defined  by  5 A.  The  kxk 
nonsingular  matrix  5  is  a  scaling  matrix  and  A  is  a  member  of  A. 
When  convenient  the  scaling  matrix  S  will  be  absorbed  into  the 
control  matrix  B2. 


What  restrictions  must  be  placed  on  the  norm  of  A  in  order  to 

ensure  that  A  +  B2AC2  is  stable?  To  answer  this  question,  define  the 
*  * 

set  e  =  e(A,  B2,  C2,  A)  and  the  subset  e  =  e(A,  B2,  C2,  A)  as, 

e(A,  B2,  C2,  A)s  {A  e  A  :  det  (si  -  A  -  B2AC2) 

=  0  for  some  s,  Res  >  0} 


e(A,  B2,  C2,  A)s  {A  e  A  :  det  (jcol  -  A  -  B2AC2) 
=  0  for  some  o),  toe  R) 
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The  set  e(A,  B,  C,  A),  which  is  similar  to  the  matricial  spectral  set 
defined  in  Reference  14,  is  the  set  of  perturbations  in  A  that  make  A 
+  B2AC2  unstable.  The  set  is  also  invariant  under  a  nonsingular 

similarity  transformation  T,  i.e.,  e(TAT  5,  TB,  CT1,  A  )  =  e(A,  B,  C,  A). 
By  analogy  with  the  frequency-domain  structured  singular  value, 
p(M,A),  we  define  the  state-space  structured  singular  value 


p(A,  B2,  C2,  A)  =  1 


supremum 
Aeg(A,B2,C2,A  ) 

0 


if  e  =  e(A,  B2,  C2,  A)  *  <(> 


if  e  =  <J) 


The  following  theorems  establish  some  useful  properties  of 
p(A,  B2,  C2  A)  and  reveal  the  connection  between  the  state-space  and 
frequency-domain  viewpoints.  The  next  three  theorems  are  proved 
in  the  Appendix.  The  first  states  that  the  supremum  always  exists  as 
a  real  number.  The  proof  of  the  second  theorem  demonstrates  that 

e(A,  B2,C2,A)  is  closed,  and  the  supremum  could  be  replaced  by  the 
maximum. 


Theorem  3.1 


If  a(A)  <  0  and  the  set  A  is  closed,  then  there  exists  a  co0  e  R 
such  that 


p(A,  B2,  Ci,  A)  <  ||  B 2  II  II C2  ||  ||  (/coo  I  -  A)"1  | 


Theorem  3.2 

If  e(A,  B2,  C2,  A)  *  0,  A  is  a  closed  set  and  a(A)  <  0,  then 

p(A,  B„  C„  A)  =  maximum  — — — 

A€^(A,B2.C J.A)  lAll 
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The  next  theorem  is  a  maximum  principle  for  p(A,  B2,  C2  A). 
Because  of  this  theorem  the  search  for  the  minimum  destabilizing 
perturbation  can  be  confined  to  the  imaginary  axis. 

Theorem  3.3 

If  £(A,  B2,  C2,  A)  is  nonempty,  A  is  a  closed  set  and  a(A)  <  0,  then 

p(A,  B,,  Gi,  A)  =  maximum  (3.3) 

Ae  eU.B^A)  II All 

The  following  theorem  establishes  the  equivalence  of  the  two 
matricial-spectral  sets,  when  M(s)  is  a  strictly  proper  (i.e.,  M(°°)  =  0) 
transfer  matrix.  It  depends  on  the  following  well-established 

determinant  identity  for  the  matrices  W  and  Z  (We  C''^,  Ze  d^) 
and  idenitiy  matrices  I(/)and  I(n)  (!(/)  e  R*x/,  I(n)  e  Rnxn)  [Reference 
15] 


det  (I  (/)  +  ZW)  =  det  (I  (n)  +  WZ) 


Theorem  3.4 

If  e(A,  B2,  C2,  A)  is  nonempty,  A  is  a  closed  set  and  a(A)<0,  then 
p(C2(/o)  I  -  Ar!B2,  A)  =  p(A,  B2,  C2,  A) 


Proof 

Since  A  is  stable,  the  resolvent  matrix  si  -  A  is  nonsingular  for 
all  s,  Res  >  0.  The  following  identity  establishes  the  claim: 

det  (sl(n)  -  A  -  B2AC2)  = 

det  (sl(n)  -  A)  det  (I(k)  -  C2(sl(n)  -  A)-1B2A)  (3.4) 

Theorem  3.4  provides  a  simple  and  direct  way  to  algebraically 
determine  the  matrices  in  Figure  1.  If  the  perturbation  matrix  SA  is 
factored  by  the  algebraic  guidelines  defined  in  this  section,  then  M(s) 
and  A  in  Figure  1  are  algebraically  determined  in  Equation  3.4. 
Moreover,  Theorem  3.4  justifies  computing  the  function  p(A,  B2,  C2, 
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A)  by  a  Monte  Carlo  eigenvalue  search  over  the  matricial  set  outlined 
in  the  paper  by  Hewer,  Kenney,  and  Klabunde  [References  8  and  9). 

The  equivalency  result  in  Theorem  3.4  indicates  that  the 
frequency-domain  and  state-space  theories  are  complimentary.  This 
adds  a  degree  of  flexibility  to  robustness  analysis  and  allows  the 

designer  to  use  whichever  approach  is  convenient.  This  is  illustrated 
in  the  following  way.  Theorems  3.3  and  3.4  provide  an  independent 
and  alternative  proof  of  the  maximum  principle  for  p  when  M(s)  is 
strictly  proper  (i.e.,  M(°°)  =  0).  Equation  3.3  is  reminiscent  of  the 
maximum  principle.  Effectively  this  principle  says  that  we  need 
only  consider  s  on  the  imaginary  axis  (i.e.,  the  boundary  of  the  right- 
half  plane)  when  computing  the  members  of  the  matricial-spectral 

set.  This  is  implicit  in  the  proof  of  the  structured-stability  theorem 
given  in  Reference  3.  A  proof,  using  the  theory  of  subharmonic 
functions,  of  the  maximum  principle  for  structured  singular  values  is 
found  in  the  paper  by  Boyd  and  Desoer  [Reference  16].  The  authors 
also  proved  that  many  basic  functions  in  control  theory  are 

subharmonic.  Theorem  3.3  is  reproved  by  using  a  simple  homotopy 
argument  that  can  be  found  in  Reference  17  and  the  equivalency 

result  Theorem  3.4. 


STABLE  REAL  PERTURBATIONS 

The  main  thrust  of  this  section  is  to  show  that  p(A,  B2,  C2  A) 
provides  a  nonconservative  measure  of  the  stability  margin  for 
stable  matrices  when  they  are  subjected  to  structured  perturbations. 
As  a  corollary  to  this  theorem,  a  decomposition  of  system  (3 . 1  )-(3 .2) 
into  a  parallel  interconnection  of  perturbed  and  unperturbed 
transfer  matrices  is  obtained.  The  subsets  A§  of  A  are  the  norm 

bounded  sets  Ag  =  (Ae  A:  1  A  |j  <8}. 

Theorem  4.1 

If  e(A,  B2,  C2  A)  is  nonempty,  A  is  a  closed  set  and  a(A)  <  0, 

then 


a(A  +  B2AC2)  <  0,  VA  e  As 

0 


1  3 


(4.1) 
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if  and  only  if 

5n(A,  B2,C2,  A)<  1 


(4.2) 


Proof 

Throughout  the  proof  let  Amax  be  a  maximal  block-diagonal 
matrix  in  e(A,  B2,  C2,  A)  whose  existence  is  guaranteed  by  Theorem 
3.3.  If  inequality  in  line  4.2  is  valid,  then  by  comparing  the 
membership  requirement  for  the  two  sets  e(A,  B2,  C2  A)  andAg,  it 
follows  that  their  intersection  is  empty.  Thus,  the  statement  in  line 
4.1  is  valid. 

Conversely,  suppose  that  4.1  is  satisfied  and  inequality  4.2  is 
false.  Thus,  it  follows  that  5>  ||  Amax  ||  and  Amax  is  a  member  of  both 
sets  e(A,  B2,  C2  A)  and  Ag,  which  means  that  4.1  is  false. 

Corollary  4.1 

If  e(A,  B2,  C2  A)  is  nonempty,  A  is  a  closed  set,  a  (A)  <  0  and 
inequality  4.2  is  valid,  then  for  all  s«  A(A  +  B2AC2)  and  s  «  A(A) 

C i  (sl(n)  -  A  -  B2AC2)_1  Bi  = 

Gu(s)+  (G12(s)A(I(k)-G22(s)Ar1  G21(s))  (4.3) 


Proof 

This  proof  relies  on  the  well-established  matrix  identity  often 
attributed  to  Sherman-Woodbury-Morrison  [Reference  7], 

(si  -  A  -  B2AC2)_1  = 

(si  -  A)-I+  (si  -  A)-1B2A(I  -  C2(sl  -  A)'‘B2A)-1C2(sI  -  A)"1 

The  block  diagram  interpretation  of  this  corollary  says  that  the 
transfer  matrix  Cj(sl  -  A  -  B2AC2)-,B1  of  the  perturbed  linear-control 
system  3. 1-3.2  with  output  feedback  matrix  B2AC2(5A  =  B2AC2)  can 
be  decomposed  into  two  parallel  systems.  In  fact  when  I  -  G22(s)A  is 
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nonsingular,  then  the  right-side  of  Equation  4.3  is  the  linear 
fractional-transformation  matrix  using  the  feedback  relation  R2  = 
AY2  in  the  two  port  T(s).  Using  the  two-port  notation,  we  obtain  the 
following  block  diagram  for  the  perturbed  and  unperturbed  system. 
Equation  4.3  and  the  block  diagram  (Figure  2)  link  the  state-space 
system  3. 1-3.2  with  the  two  port  matrix  T(s)  and  with  p(A,  B2,  C2,  A  ). 


FIGURE  2.  Canonical  Block  Diagram  for  Real  Perturbations. 


ROOT  LOCUS  AND  RELATED  EXAMPLES 


Structured  singular  value  theory  is  a  generalization  of  classical 
root  locus  and  multivariable  root  locus.  The  link  between  root  locus, 
generalized  eigenvalues,  and  structured  singular  values  is  now 
sketched.  This  link  is  used  to  compare  the  respective  stability 
margins  for  second-  and  third-order  polynomials. 

The  final  example  will  discuss  additive-error  modeling  for  the 
second-order  linear  system 


0,  1 
-coj,  -25 (0j 


C  =  [1, 


0] 


(5.1) 


with  nonnegative  damping  coefficient  %  and  positive  frequency 
variable  Wj  is  discussed.  Higher  order  transfer  functions  exemplified 
by  autopilots  [Reference  18,  p.99]  often  include  polynomials  with 
coefficients  that  are  products  of  related  variables.  For  this  reason. 
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the  additive-error  modeling  paradigm  for  5.1  is  relevant  to  many 
applied  problems. 

Any  minimal  strictly  proper  nxn  transfer  matrix  can  be 
factored  into  a  product  of  right  coprime  polynomial  matrices  N(s)  and 
D(s) 


G(s)  =  N(s)  D-1(s)  and  det  (si  -  A)=  Q1  det  D(s) 

where  Qj  is  a  nonzero  constant  [Referenced,  p.  282].  Classical  root 
locus  predicts  stability  behavior  by  computing  the  loci  in  the  complex 
plane  of  the  roots  of  the  nth -degree  polynomial  d(s)  +  A.n(s)  with 
coprime  factors  n(s)  and  d(s),  and  transfer  function  g(s)  =  n(s)d_1(s) 
as  the  scalar  parameter  X  varies. 

To  simplify  the  discussion  suppose  that  the  transfer  functions 
are  strictly  proper.  A  transfer  function  G(s)  is  exponentially  stable  if, 
and  only  if,  it  is  analytic  and  bounded  in  Res  >  0  [Reference  20].  Let 
(A,  b,  c)  be  the  minimal  realization  of  G(s)  in  controllable  conical 
form.  Using  the  basic  matrix  identity  3.4,  the  following  equation 
relates  root  locus  and  the  matricial-spectral  set  e(A,  b,  c.  A)  for  the 
triple  (A,  b,  c)  when  A  contains  only  scalar  multiples  of  the  identity 
matrix  whenever  s  e  A(A) 

det  (sl(n)  -  A  -  A.bc)  =  det  (sl(n)  -  A)  det  (1(1 )  -  A.c(sl(n)  -  A)-1b) 
=  Qj  det(d(s)- Xn(s))  (5.2) 

The  next  equation  shows  that  root  locus  is  equivalent  to  a 
generalized  eigenvalue  problem.  The  paper  by  Thompson,  Stein,  and 
Laub  [Reference  21]  discusses  multivariable  root  locus  as  a 
generalized  eigenvalue  problem.  Using  coprime  factorization  and 
letting  A  be  the  scalar  set  D(X,  X,  ...,  X),  Equation  3.4  applied  to  the 
triple  (A,  B2,  C2)  with  A  =  A.I  and  the  coprime  factors  N(s)  and  D(s) 
yields  the  following  equation  for  s  g  A(A) 

det  (sl(n)  -  A  -  XB2C2)  =  det  (sl(n)  -  A)  det  (I (k )  —  A.C2(sI(n)  -  A)_1B2) 
=  Qj  det  (D(s)-A-N(s))  (5.3) 
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Any  polynomial  subject  to  additive  perturbations  can  be 
represented  by  the  dynamical  control  system  3. 1-3.2  where  A  is  the 
companion  matrix.  The  dimension  of  B2  and  C2  and  the  location  of 
their  nonzero  coefficients  are  defined  by  and/or  will  determine  the 
perturbation  structure.  Fundamentally,  classical  root  locus  and 
multivariable  root  locus  require  a  uniform  affine  shift  of  the 
perturbed  coefficients,  while  the  structured  perturbations  only 
require  an  affine  shift.  This  observation  follows  easily  by  comparing 
Equations  3.4,  5.2,  and  5.3.  For  this  reason,  the  two  techniques  will 
generally  give  identical  robustness  measures  only  when  A  has  a  very 
special  structure.  The  following  discussion  highlights  these 
differences. 

Consider  the  third-order  stable  polynomial  p(s)  =  s3  +  a3s2  +  a2s 
+  aj  with  real  coefficients.  The  matrices 


'  0 

1 

0  ' 

O' 

A  = 

0 

.-ai 

0 

-a2 

1 

-a3. 

•b2  = 

0 

.1. 

ro  oi 


0  0 
1  1 


,c2  = 


1  0  0 
0  1  0 


define  the  root-locus  perturbation  matrix  5a  =  Xb2c2  and  the 
structured  singular-value  perturbation  matrix  5 A  =  B2SAC2  with 
scaling  matrix  S  =  -aj  ©  -a2  and  A  =  Aj  ©  A2. 


Associated  with  any  polynomial  p(s)  is  the  Hermite  matrix 
[Reference  11,  p.  465],  which  for  the  perturbed  polynomial  s3  +  a3s2  + 

a2  (1  +  A2)  s  +  aj(l  +  Aj)  is 

a  j(l  +  Aj)a2(l  +  A2),  0,  a^l+Aj) 

0,  a2(l  +  A2)a3- aj(l  +  Aj),  0 

ajO+Aj),  0,  a3 

This  matrix  can  be  used  to  test  the  stability  of  a  polynomial. 
Hermite's  theorem  says  that  a  polynomial  is  stable  if,  and  only  if, 
Hp(A)  is  positive  definite  or  equivalently,  if  its  leading  principal 
minors  are  all  positive.  Applying  this  test  leads  to  the  well-known 
condition  that  p(s)  is  stable  if,  and  only  if,  aj  >  0,  a2  >  0  a3  >  0 
and  a2a3  >  aj.  Root  locus  would  predict  a  real  margin  of  one  at  to  =  0, 


HpCA)  = 
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which  would  coincide  with  the  2-norm  of  A  =  -1  ©  -1.  The  latter 
choice  makes  the  first  principal  minor  of  Hp(A)  zero.  This  matrix  is 
clearly  a  member  of  e(A,  B2,  C2,  A),  although  it  is  not  the  minimum 
element.  The  second  leading  principal  minor  will  vanish  whenever 

the  linear  equation  Aj  =  mA2  +  m-1  is  satisfied  with  m  =  a2a3/aj. 

~  fl  -  m^ 

Since  m  >  1,  the  matrix  A  =  (0)  ©  -  is  also  a  member  of  e(A, 

v  m  y 

B2,  C2,  A)  and  has  a  2-norm  that  is  always  less  than  one.  Thus,  root- 

locus  and  the  structured  singular-value  tests  are  never  equal  even 

for  this  simple  problem. 

The  transfer  matrix  M(s)  in  Figure  1  can  be  derived  by 
computing  the  four  transfer  functions  for  all  input-output  pairs 

(Ylt  and  (Y2,  R2)  in  Figure  3. 
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By  comparison  the  state-space  approach  yields  the  algebraic 
equation  M(s)  =  C2(sl  -  A)_1B2S,  which  by  Theorem  3.4  is  the  same 


matrix,  namely  M(s)= - 

3  det(sl  -  A) 


with  A  =  Aj  ©  A2. 


The  following  discussion  shows  how  product  terms  can  be 
accommodated  by  the  state-space  formulation.  It  also  illustrates  that 
the  matrix-perturbation  structure  A  with  subset  D ,  which  is 
compatible  with  Figure  1,  does  not  include  these  product  terms.  The 
polynomial  for  system  (5.1)  is  s2  +  2£co1s  +  coj2.  The  perturbed 

polynomial  s2  +  2£coj(l  +  A j )( 1  +  A2)s  +  coj2(l  +  A|)2  correctly  models 
the  rate  of  change  and  the  product  dependency  in  the  respective 
coefficients.  The  perturbation  matrix  8A  =  B2SAC2 


B2  -= 


0  0  0 
1  1  1 


,C2  = 


0  1 
0  1 
1  0 


S  =  -  2^  ©  -  2^wi  ©  o)j 


=  A2  ©  Aj  (1  +  A2)  ©  2 , 


is  the  corresponding  perturbation  matrix  8A  in  factored  form.  Note 
that  the  product  term  precludes  the  subset  £>(Xi ,  ^2,  X.3).  However,  if 
the  set  A  included  matrix  polynomials,  then  product  terms  could  be 
accommodated.  The  stability  theorems  would  be  unaffected  but  the 
spectral  radius  lower  bound  for  p  in  Equation  2.5  would  no  longer  be 
valid.  However,  the  lower  bounds  could  be  replaced  with  the 
eigenvalues  of  matrix  polynomials.  This  example  shows  that  some 
additive-error  models  are  inherently  nonlinear,  which  complicates 
the  determination  of  M(s)  and  A  by  simple  block  diagram 
manipulation.  However,  this  example  demonstrates  that  product 
terms  can  be  directly  accommodated  by  state  space  models. 


PERTURBATION  DYNAMICS 

Both  additive  and  complex  perturbations  can  govern  additive 
errors.  Complex  perturbations  can  also  be  used  to  model  higher- 
order  unmodeled  dynamics.  Typically,  these  dynamics  have  been 
neglected  in  the  initial  design  because  they  overcomplicate  the 
design  or  because  some,  but  not  all,  of  the  general  characteristics  of 

1  9 
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the  dynamical  system  are  known.  Often  these  characteristics  are  less 
precisely  known  than  the  initial  design  dynamics. 

The  set  <7(5)  =  (G33(s):  l|G33  II  ^  <  8},  which  represents  the  pertur¬ 
bation  dynamics,  is  the  norm  bounded  family  of  exponentially  stable 

strictly  proper  transfer  matrices  ( ||  G33 II  =  supremum  ||  G33OS)  il  2X 

Res=0 


Another  theorem  in  the  Boyd-Desoer  report  shows  that  the 
latter  norm  is  subharmonic  and  thus  achieves  its  maximum  on  the 
boundary  of  s  Res  >  0.  Their  result  justifies  the  computation  of  the 
norm  only  on  the  boundary.  The  n*h -order  dynamic  system  is  a 
minimal  realization  ^3,83,03)  of  the  transfer  matrix  G33(s) 

x  =A„x„  +  B„r„  (6.1) 

33333 

y3  =C3x3  (6.2) 


The  neglected  dynamical  6. 1-6.2  system  is  interconnected  to 
3. 1-3.2  by  either  feedback  connection  defined  in  Figures  4  and  5. 


FIGURE  4.  Perturbation  FIGURE  5.  Perturbation  Feedback  With 
Feedback.  Additional  Inputs  and  Outputs. 


The  main  difference  between  these  two  systems  is  that  Figure  4 
represents  a  system  that  may  be  uncontrollable  and  unobservable, 
while  Figure  5  has  been  augmented  with  additional  inputs  and 
outputs  to  ensure  that  it  is  a  controllable  and  observable  system 
[Reference  19,  p.  446].  The  matrix  Af  of  either  interconnection  has 
the  state-variable  description 
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A,  -B2C3 

Af  = 

B3C2,  A3 

with  state  variables  x  and  x3. 

Theorem  6.1 

The  significance  of  the  following  theorem  is  that  the  structured 
singular  value  p(A,  B2,  C2,  A)  based  in  a  complex  perturbation  set  can 
be  used  to  bound  the  dynamical  family  G(8),  provided  that  the  family 
members  G33(s)  are  also  members  of  A  for  any  fixed  s  with  Res  >  0. 
The  set  identity 

e(A,  -B2,  C2,  A)  =  e(A,  B2,  -C2,  A) 

means  that  input  or  output  sign  transformations  lead  to  equivalent 
spectral  sets,  a  fact  that  is  implicit  in  the  theorem. 

Suppose  that  e(A,  -B2,  C2j  A)  is  nonempty  and  the  set  G(S)  is 
given.  If  G33(s)  e  A  for  every  s  Res  £  0,  a(A)  <  0  and  if 

8p(A, -B2,  C2,  A)  <  1  (6.3) 

then  the  nxn  matrix  I  +  G22(s)^33(s)  *s  nonsingular  for  all  s  Res  >  0 
and  any  G3 3 (s)  e  G(5).  Moreover,  the  matrix  Af  is  stable,  i.e., 
a(Af)  <  0. 

Proof 

If  there  exists  an  s0  Res0  >  0  such  that  I  +  G22  (so)  G33  (s0)  is 

singular  for  some  G33(s)  e  G(5),  it  follows  by  Equation  3.4  that 
a(A  +  (-B2)  G33  (s0)  C2)  ^  0.  The  choice  of  8  in  inequality  6.3 
implies  by  Theorem  4.1  that  a(A  +  (-B2)  AC2)  <  0  for  all  A  €  Ag. 
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But  G33(s)  is  also  a  member  of  A  for  all  s  Res  >  0.  This  contradiction 
finishes  the  first  part  of  the  proof. 

Since  A  is  stable,  the  Schur  matrix  is  defined  for  all  Res  £  0 

I,  O' 

-B3C2(sI-A)~1,  I 

so  premultiplying  s(I  ©  I)  -  Af  by  this  matrix  and  evaluating  the 
determinant  of  the  resulting  matrix  yields 


det  (s(I  ©I)  -  Af)  = 

det  (si  -  A)  det  (si  -  A3  +  B3C2  (si  -  A)"1  B2C3) 

Now  use  the  well-established  identity  in  the  Linear  System 
Model  for  Additive  Errors  section  to  obtain  the  equation 

det  (s(I  ©  I)  -  Af)  = 

det  (I  +  G33(s)  G22(s)>  (  det  (si  -  A3>  det  (si  -  A)) 

Since  A3  and  A  are  both  assumed  to  be  stable  and  I  +  G22(s)G33(s)  is 
nonsingular  all  s  Res  £0,  the  right-hand  side  never  vanishes  for  all 
Res  >0.  Thus,  the  eigenvalues  of  Af  must  occur  in  Res  <  0,  which 
completes  the  proof. 

Theorem  IV  in  Reference  20  shows  that  the  multi-input,  multi¬ 
output  system  described  by  Figure  5  is  exponentially  stable  when 
the  conditions  of  Theorem  6.1  are  satisfied.  The  methods  and 
theorems  presented  in  this  report  can  be  easily  extended  to  include 
parallel  or  series  interconnections. 


NEW  UPPER  BOUND  FOR  p 

A  new  upper  bound  for  the  smallest  destabilizing  structured- 
perturbation  matrix  8A  is  derived.  The  bound  is  theoretically 
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interesting  because  it  connects  p  with  Stewart's  [Reference  22] 
numerical  quantity  sep  (Ajj,A22). 

In  order  to  obtain  this  bound,  some  established  results  are 

reviewed.  The  separation  of  two  nxn  matrices  Aj  j  and  A2  2 
introduced  by  Stewart  is  the  quantity 

sep  (An,  A22)  =  minimum  II  AnZ-Z  A22  H 

llzi=l 

It  is  intended  to  measure  the  separation  of  the  eigenvalues  of  Aj  } 

and  A22.  In  fact  when  Aj  j  and  A22  are  normal  matrices,  then 

sep(A]j,  A22)  is  the  minimum  distance  between  the  eigenvalues  of 

A  j  j  and  A22  for  the  2-norm.  An  important  and  useful  property  of 
the  sepCAj],  A22)  is  its  relation  to  the  Sylvester  equation 

AuX-XA22  =  W  (7.1) 

A  central  fact  in  dealing  with  (7.1)  is  that  it  is  equivalent  to  a  linear 
system  of  order  n2  [Reference  11,  p.  414]  Px  =  co  where  P  =  I  <8>  An  - 
A22T  ®  I  and  x  =  Vec  (X),  w  =  Vec  (W).  Here  ®  indicates  Kronecker 
product  and  Vec  (E)  is  the  n2  square  vector  formed  by  stacking  the 
columns  of  the  nxn  matrix  E.  For  the  Frobenius  norm  the  following 
equation  has  been  known  for  some  time  [References  23  and  24] 

0mir(P)  -  „  ™in  tl  Px  II  2  =  sepF  (A  1 1 ,  A22)  (7.2) 

14  =  1 

where  amin(-)  denotes  the  minimum  singular  value.  Recently,  Hewer 
and  Kenney  [Reference  25]  and  Kenney  and  Hewer  [Reference  26] 
established  the  new  bound  for  stable  matrices  A 

sep2(AT,  -A)=  IIHH21  (7.3) 


with 


AtH  +  HA  =  -I 


(7.4) 
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Thus,  for  both  norms  the  separation  between  two  matrices  can 
be  easily  computed.  Moreover,  the  two  expressions  7.2  and  7.3 
provide  a  new  link  via  the  following  theorem  between  Lyapunov 
equations  and  structured  singular  values. 

Theorem  7.1 


If  a(A)  <  0  and  A  is  closed,  then  p(A,  B2,  C2,  A  )  < 


2  H  B2 II  1 C21I 

sep  (At,  -A) 


Proof 

For  any  A  €  e(A,  B2,  C2,  A)  the  matrix  A  +  B2  A  C2  has  an 
eigenvalue  on  the  imaginary  axis.  Hence,  the  equation 

(A  +  B2AC2)TX  +  X  (A  +  B2AC2)  =  0 

has  a  nonzero  matrix  solution  X.  Rewriting  we  have 

AT  X  +  X  A  =  ((B2AC2)TX  -  X  B2AC2) 


thus 


sep  (A,  -A)  < 


AtX  +  X  A 


<  2 II  B2AC2||  <  2  B2  C2 


This  gives 


1  ^  2  11  B2  H  11  C;  H 

II A  ||  sep  (AT, -A) 

Then  taking  the  maximum  over  all  A  e  e(A,  B2,  C2,  A)  gives  the  result. 


DISCRETE-TIME  CASE 

All  of  the  above  results  are  stated  for  the  continuous-time  case. 
All  of  the  proofs  that  are  purely  algebraic  and  depend  on  simple 
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properties  of  rational  matrices,  determinants,  and  matrices  apply 
equally  well  to  the  discrete-time  case  with  some  modification.  The 
Laplace  transform  is  replaced  by  the  Z-transform.  The  role  of  the 
left-half  imaginary  axis  and  right-half  plane  is  taken  by  the  interior 
boundary  and  exterior  of  the  unit  disk,  respectively.  The  stability 
condition  a(A)  <  0  is  replaced  by  the  condition  p(A)  <1,  and  the 
stability  tests  are  rewritten  with  ei10  replacing  jo).  Moreover,  the 
arguments  in  the  continuous  case  can  be  replaced  with  the 
corresponding  proofs  and  bounds  for  the  discrete  case  that  are  found 
in  Reference  17.  The  results  of  the  previous  section  can  be  derived 
using  the  Cayley  matrix-operator  mapping  (Reference  27)  that  maps 
the  stable  discrete  Lyapunov  equations  onto  the  stable  Lyapunov 
equation  (7.4). 

DE  GASTON  EXAMPLE 

The  example  in  Figure  6  is  found  in  the  thesis  of  de  Gaston 
[Reference  4],  and  captures  some  of  the  essential  features  of  an 
autopilot  [Reference  18].  This  example  illustrates  Doyle’s  claim  that  a 
linear  system  can  be  rearranged  to  match  Figure  1  and  demonstrates 
state-space  modeling.  The  nominal  parameter  values  are  K  =  800, 
z  =  2,  p  =  10,  co2  =  6,  o>3  =  4,  0)3  =  4,  0)2  =  6. 


FIGURE  6.  Compensator  and  Dynamical  System. 

The  block  diagram  in  Figure  6  was  converted  to  a  block  diagram  with 
scaled  perturbations  (Figure  7).  The  scaling  was  selected  so  that  a 
value  of  ±  1  for  the  perturbation  represents  100  percent  of  the 
nominal  parameter. 
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FIGURE  7.  Compensator  and  Dynamical  System  With  Output  Scaling. 

Aj  (i  =  1,  2,  3)  represent  real  perturbations,  and  toj  (i  =  2,  3)  are  the 
nominal  poles  and  scaling  parameters  such  that  the  perturbation 
represents  a  percentage  of  the  nominal  parameter. 

Figure  6  is  transformed  to  the  desired  form  by  rearranging  the 
loops  at  the  inputs  to  the  perturbation  blocks  (Figure  7).  Treating 
the  inputs  and  outputs  of  these  blocks  as  separate  system  inputs  and 
outputs,  respectively,  one  can  then  use  conventional  block  algebra 
and  Figure  8  to  obtain  the  desired  M  matrix  by  evaluating  all  nine 
transfer  functions  in  Figure  9.  The  original  system  inputs  and 
outputs  are  no  longer  represented  in  Figure  8. 


FIGURE  8.  Open-Loop  Perturbation. 
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y3 

Ri 

r2 

R3 

FIGURE  9.  Input-Output  Matrix. 

Figure  10  is  the  M  matrix,  which  is  found  in  References  4  and  5 
with  o>2  =  o>3=  1,  for  the  system  defined  in  Figure  6. 


R  1 


«2 


r3 


-K(s  +  z) 

K  (s  +  z) 

Kfe  +  z)  (s  +  o)2) 

d 

d 

d 

s(s  +  p)  (s  +  0^)  0^ 

-s(s  +  p)  (s  +  0^)  0^ 

K  (s  +  z)  0^ 

d 

d 

d 

s(s  +  p  )  0)3 

— s(s  +  p)  0)3 

s(s  +  p)  (s  +  0)2)0)3 

d 

d 

d 

d  =  s  (s  +  0)2)  (s  +  0)3)  (s  +  p)  +  K(s  +  z) 
FIGURE  10.  Sturctured  Feedback  Matrix. 
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The  controllable  and  observable  linear-dynamical  system  for  Figure 
6  can  be  derived  with  the  aid  of  the  state  variables  xj,  X2,  X3,  X4  that 
are  identified  in  Figure  11. 


FIGURE  11.  State  Variables. 


O 

X 


1  =  -pxj  +  (z  -  p)e 


yi  =  xj  +  e 


*2  =  Yi 

X3  =  — (©2  +  80)2)  X3  +  (K  +  8K)  Xj 

X4  =  -(<03  +  80)3)  X4  +  X3 


e=u-y  =  u-x4 

These  equations  yield  the  following  dynamical  system  (A,  Bi,  Ci) 


A  = 


B,= 


-p  0  0 

p-z" 

1  0  0 

-1 

§ 

1 

u 

0 

0 

0  0  1 

-G)j, 

1  matrices 

/ - 

0  0 

0  0 

0  0 

B,  = 


_\ 

z-p 

1 

0 

v  0  ; 


1  -1  0 

VO  0  -1 


C2  = 


"0  1  0  0" 
0  0  10 
,0  0  0  1, 


Cl  =  (0,  0,  0,  1) 


=  Ai  ©  A2  ©  A3 
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and  scaling  matrix  S  =  K  ©  (02  ©  0)3. 

The  perturbation  matrix  5A  with  8K  =  KAi,  8t02  =  -g)2A2,8co3  =  -CD3A3 

r0  0  0  0  ' 

0  0  0  0 

8A  = 

0  KAj  — 002A2  0 

^0  0  0  -cogAj, 

factors  into  the  product  8A  =  B2SAC2. 

In  order  to  determine  the  state-space  equivalence  of  M(s)  in 
Figure  10,  the  inverse  of  the  matrix  (si  -  A)  must  be  determined. 
Let  (si  A)ij  denote  the  (n  -  1)  x  (n  -  1)  matrix  formed  by  deleting 
row  i  and  column  j  from  A.  Define  the  cofactor  matrix 

cof(sI  -  A)  =  (-l)i+J  det  (si  -  A)ij,  i,j  =  1,  ...n 

then  si  -  A  has  an  inverse  equal  to  det1  (si  -  A),  times  the  transpose 
of  the  matrix  obtained  from  (si  -  A)  by  replacing  each  element  by  its 
cofactor  whenever  s  t.  A(A).  For  the  de  Gaston  example,  it  is  easy  to 
verify  that 


d  =  det  (si  -  A) 

where  d  is  defined  by  the  equation  in  Figure  10. 

The  transpose  of  the  cofactor  matrix  for  (si  -  A)  is 

'K+s(s+C02)(s+G>3),  -KCz-p),  -s(z-p),  s(s+C02)(p-z)' 

(s+G^Xs+Cflj),  (s+pXs+tojXs+cog),  -(s+z),  -(s+o^Xs+z) 

K(s+co3),  K(s+pXs+co3),  sfs+pXs+wj),  -K(s+z) 

,  K,  KCs+p),  s(s+p),  s(s+pXs+a2), 

The  matrix  M(s)  =  SC2(sI  -  A)’,B2  is  identical  to  the  transfer 
matrix  in  Figure  10.  The  transfer  matrix  with  input  scaling  Mi(s)  = 
C2(sl  -  A)_1B2S  is  not  equal  to  M(s)  with  output  scaling.  However,  by 
the  fundamental  determinant  identity,  the  polynomials  in  the  three 


29 


NWC  TP  6971 


variables  Aj,  A2,  and  A3,  defined  by  the  equations  det  (I  -  Mi(s)A) 
and  det  (I  -  M(s)A),  are  identical.  This  observation  leads  to  the 
following  theorem  about  input-output  scaling. 

Theorem  8.1 

If  a  (A)  <  0  and  e  (A,  B2,  C2,  A)  are  nonempty,  then  for  any 
nonsingular  transformation  T  that  commutes  with  the  members  of  A 

P  (A,  B2T,  C2,  A)  =  p  (A,  B2,  TC2,  A) 


CONCLUSION 

A  state-space  theory  of  structured  uncertainty  that  offers  a 
complementary  and  alternative  perspective  on  Doyle’s  structured- 
value  analysis  was  presented.  This  theory  brings  a  degree  of 
flexibility  to  robustness  modeling  that  allows  a  designer  to  use  either 
a  frequency-domain  or  a  state-space  viewpoint.  A  second-order 
oscillator  model  was  used  to  illustrate  the  efficiency  of  state-space 
uncertainty  modeling  when  product  terms  are  present.  The  theory 
was  compared  with  root  locus  using  a  third-order  polynomial  model. 
Although  not  developed  in  this  report,  eigenvalue  formulas  that  can 
be  combined  with  a  Monte  Carlo  search  strategy,  which  offers  an 
alternative  means  for  the  computation  of  p,  were  presented.  Finally, 
the  state-space  perspective  offered  an  attractive  alternative  tool  in 
uncertainty  modeling  and  p  computations.  A  new  stability  theorem, 
characterizing  a  family  of  neglected  dynamics  that  can  be  safely 
connected  to  the  original  system  by  feedback  without  destabilizing  it, 
was  provided.  The  family  is  confined  to  a  ball  whose  radius  is 
determined  by  a  structured  singular  value  inequality. 


ADDENDUM 

After  this  report  was  written,  we  received  a  copy  of  "Robust 
Control  of  Multivariable  and  Large  Scale  Systems,"  by  Andrew 
Packard  and  John  C.  Doyle,  Honeywell  Systems  and  Research  Center, 
Minneapolis,  Minn.,  23  March  1988,  Report  No.  F49620-86-C-001  . 
Their  report  gives  a  fairly  complete  introduction  to  the  structured 
singular  value  (p),  and  details  some  robustness  theorems  for  linear, 
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time-varying,  and  nonlinear  systems.  The  Packard-Doyle  report 
identifies  a  state-space  model  for  structured  singular  values  similar 
to  the  input-output  viewpoint  introduced  in  this  report.  Especially 
noteworthy  is  the  Packard-Doyle  description  of  how  parametric 
uncertainty  in  state-space  models  can  be  rearranged  into  the 
framework.  Their  state-space  models  and  requirements  on  the  set 
are  identical  to  the  hypothesis  advanced  in  Theorem  6.1.  Their 
theorem  subsumes  the  latter  results  and  is  an  important  and  useful 
contribution  to  the  higher  order  unmodeled  dynamics  perturbation 
literature. 

A  brief  summary  of  some  of  the  results  in  this  report  are  in  the 
publication  "A  State  Space  Model  for  Structured  Singular  Values," 
27th  IEEE  Conference  on  Decision  and  Control,  Austin,  Tex.,  7-9 
December  1988,  pp.  2144-2147. 
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Appendix  A 

PROOF  OF  THEOREMS  3.1,  3.2,  AND  3.3 

This  appendix  proves  that  the  first  three  theorems  in  the 
section  of  this  report  entitled  Linear  System  Model  for  Additive 
Errors  are  sketched. 

Proof  of  Theorem  3.1 


Since  A  is  stable  the  resolvent  matrix  jcol  -  A  is  nonsingular. 

Thus,  it  follows  from  Martin  and  Hewer  [Reference  17]  that  there 

* 

exists  a  (dq  e  R  such  that  for  every  Ae  z  (A,  B2,  C2,  A) 


I|b2AC2||  >  ||  (jojQ  -  A)-1 1|_1 


Let  P  be  the  following  set  of  real  numbers 


P  = 


:  Ae  e(A,  By  C2,  A)l 


(A.l) 


(A. 2) 


By  Equation  A.l  the  set  P  is  bounded  above;  so  by  a  fundamental 
axiom  [Reference  28,  p.  7]  of  the  real  number  system,  P  has  a 
supremum.  These  two  results  clearly  establish  the  claim. 

Proof  of  Theorem  3.2 

The  set  P  defined  by  Equation  A. 2  is  bounded  above  and  below, 

and  thus  has  an  infimum  and  supremum  [Reference  28,  p.  7].  If  the 

set  P  has  a  finite  number  of  points,  then  the  claim  is  obvious. 
Otherwise  P  is  a  bounded  infinite  subset  of  the  positive  real 
numbers.  It  follows  from  standard  arguments  in  Apostol  [Reference 
28,  pp.  64-73],  that  the  set  P  is  closed  on  the  right  because 
Ostrowski's  theorem  [Reference  15]  on  the  continuity  of  eigenvalues 

can  be  used  to  show  that  limits  of  convergent  sequences  in  e(A,  B2, 

C2,  A)  must  belong  to  e(A,  B2,  C2,  A)  since  A  is  a  closed  set. 
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Proof  of  Theorem  3.3 

Proof  of  Theorem  3.3  is  found  in  Martin  and  Hewer  [Reference 
17]  and  is  included  for  completeness. 

Since  a(A  +  B2AC2)  ^  0  for  some  Aee(A,  B2,  C2,  A),  there  exists  X 
e  C  with  Re  X  £  0  such  that  XI  -  A  +  B2AC2  is  singular.  If  X  =  j  to,  for 
some  to  e  R,  let  8  =  1.  Otherwise,  define  W(z,8)  =  zl  -  A  -  5  B2AC2  (z  <= 
C,  8  e  [0,  1]),  and  observe  that  the  determinant  of  W  is  a  monic 
polynomial  in  z  of  degree  n  whose  coefficients  depend  continuously 
on  8  e  [0,  1].  Thus,  as  8  varies  between  0  and  1,  the  roots  of  det 
W(z,8)  will  trace  continuous  curves  in  the  complex  plane.  One  of 
these  curves,  say  (3:[0,  1  ] — >C,  will  have  X  as  its  endpoint  when  8=1. 
So  we  have  (3(0)  e  A(M),  Re((3(0))  <  0,  and  (3(1)  =  X.  Since  the  interval 
[0,  1]  is  connected  and  (3  is  continuous,  the  image  (3  e  [0,  1]  is  also 
connected,  and,  therefore,  intersects  the  imaginary  axis  at  some  point 
jco.  Consequently,  there  must  be  a  value  of  8  e  [0,  1]  call  it  80,  such 
that  P(8o)=yco.  Thus,  0  =  det  W((3(8o),80)  =  detW(jco,SQ),  which 
means  ycol  -  A  -  8oB2AC2  is  singular.  Since  the  norm  obeys  the 
inequality,  8Q  || B2 AC2  ||  <  || B2AC2  | ,  the  maximum  will  occur 

somewhere  on  the  boundary  of  the  imaginary  axis  in  the  complex 
plane. 
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